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ABSTRACT 

We present an analytic calculation of the thermonuclear depletion of lithium in 
contracting, fully convective, pre-main sequence stars of mass M <^ O.SMq. Previous 
numerical work relies on still-uncertain physics (atmospheric opacities and convection, 
in particular) to calculate the effective temperature as a unique function of stellar mass. 
We assume that the star's effective temperature, Teg, is fixed during Hayashi contraction 
and allow its actual value to be a free parameter constrained by observation. Using this 
approximation, we compute lithium burning analytically and explore the dependence 
of lithium depletion on Teg, M, and composition. Our calculations yield the radius, 
age, and luminosity of a pre-main sequence star as a function of lithium depletion. This 
allows for more direct comparisons to observations of lithium depleted stars. Our results 
agree with those numerical calculations that explicitly determine stellar structure during 
Hayashi contraction. In agreement with Basri, Marcy, and Graham (1996), we show 
that the absence of lithium in the Pleiades star HHJ 3 implies that it is older than 100 
Myr. We also suggest a generalized method for dating galactic clusters younger than 
100 Myr (i.e., those with contracting stars of M O.OSMq) and for constraining the 
masses of lithium depleted stars. 

Subject headings: open clusters and associations — stars: abundances — stars: evo- 
lution — stars: fundamental parameters — stars: low-mass, brown dwarfs — stars: 
pre-main sequence 
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1. Introduction 

Lithium depletion in gravitationally contracting, fully convective stars of mass M < Mq 
(Hayashi & Nakano 1963; Bodenheimer 1965) and M = Mq (Weymann & Moore 1963; Ezer 
&; Cameron 1963) has been studied for over thirty years. Motivated by lithium observations of 
main sequence stars in young clusters and the halo (e.g., Soderblom 1995), many have calculated 
lithium depletion for stars of M ;^ O.5M0 both before and during the main sequence (Vauclair 
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et al. 1978; D'Antona & Mazzitelli 1984; Proffitt & Michaud 1989; Vandenberg & Poll 1989; 
Swenson, Stringfellow, &; Faulkner 1990; Deliyannis, Demarque, & Kawaler 1990). For stars with 
M ^ 0.5Mq, most lithium burns after convection has halted in the stellar core. Accurate predictions 
of lithium depletion then depend on the temperature at the bottom of the retreating convective 
zone. The location of the convective/radiative boundary and the amount of mixing across it depend 
on opacity, treatment of convection, and rotation; proper handling of these effects remains an open 
question. 

Lower mass (M ^ 0.5Mq) stars are fully convective during lithium burning, which occurs 
before the star reaches the main sequence.^ The effective temperature T^s of a fully convective star 
determines the contraction rate and is found by matching the entropy at the interior to that at the 
photosphere (see Stahler 1988 for a review of Hayashi contraction). For these low mass stars with 
Teff ^ 4000 K, the opacities are still uncertain, and there are still debates about the treatment of 
convection. These uncertainties have motivated many numerical calculations of lithium depletion 
(Pozio 1991; Nelson, Rappaport, & Chiang 1993; D'Antona & Mazzitehi 1994 (hereafter DM94); 
Chabrier, Baraffe, & Plez 1996). The differing input physics results in different Tpfj's for the 
same stellar mass. For example, DM94 found that, depending on the opacities used, T^g ranges 
from 3350 K to 3640 K for a O.2M0 star. Most calculations agree that Teg remains approximately 
constant during the fully convective contraction phase. 

In this paper, we present a different approach to calculating pre-main sequence (or pre-brown 
dwarf) lithium depletion in gravitationally contracting, low mass (M ^ 0.5Mq) stars. Rather than 
calculate Teg, we calculate the dependence of lithium depletion on Tefj. Given Tgfj, M, and mean 
molecular weight /i, we can reproduce the results of prior works. Moreover, our approach allows 
the inferred effective temperature to be used directly in analyzing lithium depletion observations. 
Efficient convection throughout the star allows it to be modeled as a fully mixed n = 3/2 polytrope. 
Our analytic calculations then yield the age, radius, and luminosity at a given level of lithium 
depletion. Our results (equations |11] and II^) can be used to constrain the mass and age of a star 



from its lithium abundance. We apply these same techniques to depletion of beryllium and boron 
in Ushomirsky et al. (1997). 



2. Contraction and Lithium Burning in Fully Convective Pre-Main Sequence Stars 

A number of authors have undertaken detailed numerical evolutionary calculations of con- 
tracting pre-main sequence stars (see Burrows &: Liebert 1993 for a review). Low mass stars 
(M ^ 0.5Mq) remain completely convective at least through the end of lithium burning (DM94). 
For an ideal gas (P = pNAksT/ fi where Na_ is Avogadro's number) the adiabatic relation is 



The interstellar lithium abundance is so low {NLi/NH ~ 10"'') that the energy released by its fusion does not 
appreciably slow stellar contraction. Throughout this paper, we only consider ^Li, as ®Li always depletes first and is 
less abundant in the local ISM (Lemoine, Ferlet, & Vidal-Madjar 1995). 
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P oc p^^^, so that an n = 3/2 polytrope describes the stellar structure. During later stages of 
contraction, electron degeneracy modifies the equation of state by reducing the central tempera- 
ture below its non-degenerate value. However, partial degeneracy does not alter the polytropic 
structure. The central density and temperature are 



\MqJ \ R J ^ ' ■ \0.6J\MqJ\R 

where 

_ pNAhsT 

Pes = p < P (2) 

accounts for the deviation in electron pressure from that of an ideal gas. We neglect Coulomb and 
ionization corrections to the equation of state. Except during deuterium burning, gravitational 
contraction powers the star's luminosity, L = AirR'^asBTefi^ = —{3GM'^/7R'^){dR/dt), which is 
independent of the degree of degeneracy. 

We define the time coordinate t under the assumption that Tefj is constant during contraction 
from a formally infinite radius. This time therefore differs from chronological age because of the 
deuterium burning phase and the initial radius on the theoretical stellar birthline ( ptahler 1988 ). 



Lithium depletion, however, occurs long (10-100 Myr) after these events, so that if the effective 
temperature used is correct at the time of lithium depletion, t differs only slightly from chronological 
age. Gravitational contraction then gives the stellar radius and luminosity as functions of time as 
R f M \2/3 / 3000 KX^'/V Myr 

_L = 5.25xlO-f^V'V^V'VMy£)'", (4) 

Lq \0.1MqJ V3000Ky \ t J ^ ^ 

We define the contraction timescale in terms of the central temperature, 
R /3000K\^ /O.1M0\ /0.6\ V 

which we then compare to the timescale for lithium destruction due to the reaction ''Li(p, a) ^He. 
Over the range of temperatures (Tg = (T/10^ K) < 6) appropriate for this work, the reaction rate 
is Na{(tv) = SfscrT^^^^exp{-aTQ^^^)cm'^s-^g~^ where 5 = 6.4 x 10^° and a = 84.72 (Caughlan 
& Fowler 1988) and fscr is the screening correction factor (Salpeter & Van Horn 1969). Raimann 



(1993) recently discussed new experimental results at low energies (~ 11-13 keV) and updated S 
to 7.2 X 10^°. 

Lithium is depleted when the local nuclear destruction time, t^cst = mp/X p{av) {X is the 
hydrogen mass fraction and rup is the proton mass), becomes comparable to tcont- Although the 
full calculation involves integrating over the star (§3), one obtains a flavor of the calculation by 
evaluating tjest at the center of the star (with /scr = 1 and p^g = p) and equating it with tcont- 
This gives a relation for the central temperature Tcq = Tc/W^ K at the time of depletion, 

32.9 + ln(5) - 31n f - 41n f - 61n f-^) + ^ InT,,. (6) 



Tce^/^ ' ^ ' \0.1MqJ V3000Ky V0.6y 3 



4 



Bildsten et al. 



For example, if Tcfj = 3500 K and M = 0.5Mq, then ^Li depletes when Tcq = 3.04. The temperature 
is only slightly different for other masses. For these temperatures, the reaction rate is extremely 
temperature sensitive (oc T'^^), so that t^est decreases rapidly as the temperature increases. Hence 
the transition from no depletion to full depletion occurs rapidly as long as the central temperature 
reaches these values. This does not occur for very low masses O.O6M0), because electron 
degeneracy pressure dominates the equation of state (Pozio 1991; Magazzii, Martin, & Rebolo 
1993; Nelson et al. 1993). Also, stars with M ^ O.SMq develop radiative cores prior to depletion. 
We only discuss stars safely between these limits. Since this depletion occurs on a contraction 
timescale (10-100 Myr), the central temperature during depletion is higher (~ (3-4) x 10^ K) than 
the temperature needed (2.4 x 10^ K) at the base of a main sequence star's convective zone to 
deplete ^Li while the star is on the main sequence. 

The narrow range of burning temperatures allows the solution of equation ^ to be approxi- 
mated as a power law. For the non-degenerate case, Tc is oc M^/^Tefj^/^//"^/^ at the time of depletion, 
and tdcpi is oc Tcff^'/^Af"^/®. The radius at depletion, R oc M'^/^ fi^/'^T^fi'^/^ , is relatively insen- 
sitive to Tcfj but nearly proportional to mass. These scalings, which are modified only slightly in 
the full calculation (§3), provide an intuitive picture of lithium burning and a means to evaluate 
observations (§4) once the prefactors are known. 



3. Full Calculation of Lithium Depletion 

Since lithium burns near the center of the star, depleting lithium throughout the star requires 
mixing lithium-poor fluid outward and lithium-rich fluid inward. For efficient convection, the 
mixing timescale (~ 10-100 yr) is much shorter than both tcont and tdest! convective mixing keeps 
the lithium to hydrogen ratio, /, fixed throughout the star as the total lithium content decreases. 
Thus we write the global depletion rate as 

df Xfr^ 

= / p{av)dM. (7) 

dt rup Jo ' ^ ^ 

Changing to spatial variables and using the thermonuclear rate defined earlier, we obtain 

d . „ 47rX 



dt 



In/ = --i^ / rV'5/s„T6-2/3exp ( dr. (8) 



The temperature sensitivity of the nuclear reaction allows us to expand T and p about their central 
values. Integrating the lowest order terms (see Ushomirsky et al. 1997 for a discussion of the small 
errors introduced by this approximation) yields 



where a = oTcq^^^^ is a convenient representation of the central temperature. Using equation (|5|) 
and the fact that da/dR = a/3R, we characterize the stellar state by a and integrate from a = 00 
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and initial abundance /o to find the depletion W as a function of a, 

where g{a) = a-^VSg-a _ 29r (-37/2, a) and r(-37/2, a) is an incomplete gamma function. This 
transcendental equation is similar to that found from the simple arguments in §2. For a given 
depletion W, the solution a(VF) determines Tc. The radius, luminosity, and age then follow from 



equations (|3|)-(|5D. Our treatment of degeneracy and screening in equation (10) is approximate in 
the sense that we neglect the slight variation of /igfr and /scr with a. Instead, we evaluate their 
values at the time of depletion to solve the resulting transcendental equation. We discuss this 
process in more detail in a forthcoming paper (Ushomirsky et al. 1997). 

Although solving the depletion equation (10) numerically is straightforward, an analytic fitting 
formula for the time (and therefore the radius and luminosity from equations Q and Q) as a 
function of depletion is also helpful. The age at a given depletion (we use Raimann's [1993] rates; 
using Caughlan's and Fowler's [1988] rates changes tdepi by no more than a few percent) for 2000 K < 
Teff < 4000 K, O.O75M0 <M < O.SM©, and 0.65 <X < 0.75 is 

, /o.iM0\°-^iV3oooK\3-5i /o.ex^-^V ^ / X 

tdepi = 54.1 —— — i^) (11) 



1 + 0.117 



M J \ Tcs J \ H J Vln2 

O.lM0\6-39/ Teff \ 0-828 /Qgxll.S /^x 0.204- 



M J VSOOOKy \ fi J Vln2 



Myr. 



The first line of equation ( |ll|) follows from the simple arguments of §2, while the expression in 
square brackets accounts for the onset of degeneracy in lower- mass stars. A simpler fit, 

/O.lM0\O-663 /3000K\3-5V0.6\2°9 / ^ , x 
,,„ = 50.7(^) (— ) (-) (^) Myr, (12) 

is obtained if degeneracy is not important during depletion (O.2M0 ^ M ^ 0.5Mq). These fits 
agree with Chabrier et al. (1996, Table 1) for masses O.O8M0 < M < O.5M0 at 50% (99%) depletion 
to better than 4% (8%), 7% (17%), and 12% (25%) for the radius, luminosity, and age, respectively.^] 
We also reproduce the results in Tables 5-8 of DM94 for O.O7M0 < M < 0.3Mq to better than 
3%, 7%, and 15% for the radius, luminosity, and age, respectively. 

The strong temperature dependence of the burning rate makes T^. at depletion very insensitive 
to our approximations. The central temperature is thus the most accurately determined quantity. 
Indeed, our Tc's at a given depletion level deviate from those reported in Tables 5-8 of DM94 by 
less than 2% (Chabrier et al. [1996] did not report Tc). Our neglect of Coulomb and ionization 
effects ( Ushomirsky et al. 1997] ) introduces a small error in mapping Tc to R (roughly, Tc oc M/R). 



^Percentages reported here are maximum deviations of our results (equations and [Q) from the stated 
references. 
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These errors then propagate into the age and luminosity determination at a given depletion level. 
Since the age is oc (see equation ) , the relative deviations in t are roughly three times larger 
than those in R. Moreover, any changes in Teg during contraction cause additional uncertainty 
(typically up to 5%, see Ushomirsky et al. 1997) in calculating how long a star takes to contract to 
a certain radius. 

In summary, the absolute error of age determination in our calculation is no more than 20-25% 
under any conditions, and the average discrepancy with published results is about half that amount. 
In light of these uncertainties, it is important to note that, because the age at a level of depletion 
is t oc the typical 5% observational errors in Teg lead to uncertainties in age that are 

comparable to the deviations of our results from the detailed numerical ones. Our results therefore 
are adequate for observational work, and have the added advantage of enabling the observer to use 
the derived Teg and L in conjunction with lithium observations to constrain the stellar mass. This 
is in contrast to using detailed evolutionary tracks to infer the mass from the star's location on the 
HR diagram. As we show in the following section, the inferred mass can be in conflict with lithium 
depletion. In addition (insofar as Teg is nearly constant during contraction), our results can serve 
as benchmarks in comparing and evaluating detailed numerical calculations. 



4. Comparison to Observations of Lithium Depleted Pre-Main Sequence Stars 

The radii of pre- main sequence stars are typically inferred from T^s and L, but (except in 
binaries) masses must be estimated by relating colors to a grid of computed evolutionary tracks. 
There are potentially significant uncertainties, both observational and theoretical, associated with 
inferring masses and radii of stars from these observations. On the other hand, the presence or 
absence of lithium is an indirect indicator of the central temperature of the star. Since lithium does 
not deplete until a characteristic central temperature (Tc ~ 3-4 x 10^ K) is reached, the lithium 
abundance can be used to constrain the mass to radius ratio (cf. equations ^ and [|^). Age 
determination is, in some sense, secondary to measuring the radius. 

Since the central temperature is cx M/R, the additional constraint provided by the lithium 
abundance measurements (that are sensitive to the derived Teff) is most easily exploited by plotting 
stellar radii versus masses. Figure |l| is a compilation of published mass and radius values of pre- 
main sequence stars where lithium abundances have been estimated or constrained. We divide the 
stars into high (small circles), low (large circles), and undetectable (triangle) lithium content bins 
(see caption). Also plotted are lines of constant depletion, which correspond to nearly constant 
central temperature, i.e., R oc M. These lines are relatively insensitive to the effective temperature 
and divide the graph into an undepleted, a depleting, and a depleted region. The narrowness of the 
depleting region results from the reaction's strong temperature sensitivity. We do not continue our 
lines to high masses where the stars develop radiative cores prior to depletion, nor to low masses 
where the stars become strongly degenerate. 
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A compilation of published values ([Burrows. Hubbard, fc Lunine 198£ : Magazzii, Martin, & Rebolo 
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sequence stars. Small circles represent stars reported to have [Lij = L'Z + log(I\Li/l\u) > 1.82, which corresponds to an 
abundance of more than 5% of the Population 1 value ([Li] — 3.1). Large circles correspond to stars with [Li] < 1.82. 
Triangles, both filled and empty, indicate stars with no firm lithium detection: either no lithium line is detected, or 
only an upper bound is reported for A''Li. Where applicable, we have indicated repor l^ed uncertainties in mass and 
radius. The two solid lines are contours of constant depletion calculated from equation (111) with fi — 0.6 and X — 0.7 
and with an effective temperature law TeS = AOOO{M / Mq)^^"^ K; the position and shape of the depletion contours are 
insensitive to these chnices . The solid lines end at M = O.SM©; above this mass stellar cores are not convective at 
the depletion time ( DM94 ). At lower left are a series of points (filled right-pointing triangles) for which no lithium 
was detected and no mass was reported (cf. Table ^). Taking the non-detection of lithium to mean that these stars 
are more than 95% depleted, we solv e for th eir minimum masses. The dotted line is the theoretical zero-age main 
sequence for O.O8M0 <M< 0.5Mq ([DM94|). 
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A few stars appear depleted but lie above the depletion lines. Since central burning cannot 
have caused such a low lithium abundance, either the star was born with an anomalously low 
abundance or the reported masses and/or radii are in error. A challenge to our assumption of 
complete mixing would be a detection of lithium in a star lying beneath the depletion lines. We 
have found no such measurement for M < 0.5Mq. We also display in Figure [l| those stars for which 
no lithium is detected and no mass is reported (filled triangles). We interpret the lack of lithium 
to imply that more than 95% has burned and use the inferred radii to constrain their masses (see 
Table 1). Determinations of Tgff from the observed colors are still highly uncertain and depend on 
the temperature scale and model atmosphere used. We have not attempted to convert observed 
colors to a single temperature scale, but have just used the quoted results. However, our figure 
depends only on the value of Tgfj and not the physics needed to infer it. 

Figure 1 also clearly shows that there is a dearth of observations of pre-main sequence stars 
near the depleting region with masses O.IMq < M < 0.5Mq. Observations in this mass range are 
important, as they would test the mixing assumption and may yield, independently of atmospheric 
physics, Tgfj as a function of mass. 

As mentioned above, there are significant uncertainties in relating the observed colors to ef- 
fective temperatures for low-mass stars, and in using the inferred Tgfj's and evolutionary models to 
determine ages of stars. Nevertheless, we can still bound the stellar age and mass without relying 
on exact knowledge of T^g, as long as L is well determined and lithium is unobserved^ down to 
some detection limit Wo- For a trial value for the age of the star, t, the contraction equation (Q) 
yields M oc t^/^/Tgfj. Substituting this relation into the depletion formula, equations (11) or (|l2[). 



we obtain the time t^epi at which this star is depleted to Wo- For consistency with the lack of 
lithium, we demand t > tdcpb with equality yielding the lower limit on age for a given Tpg. For 
high-mass stars, t is oc T^'^ , and the upper bound on Tefj determines the minimum age. For lower 
mass stars, t(Tcff) always has a minimum (due to the effects of degeneracy), 



,„, = 58,3(^— (-) (idj <1^> 

which is the lower bound on age independent of Tefj for —3.50 < log(L/L0) < —2.32. Because the 
minimum in tdepi is due to the effects of degeneracy, this technique is independent of Teg only for 
low-mass stars; this is fortunate because for these stars Tcfj is more uncertain. 

We apply equation (|l3|) to the lithium-depleted star HHJ 3 in the Pleiades (log[L/Lo] = 
—2.78 lb 0.05 [ Basri et al. 199^ ]). If the absence of ''Li in its spectrum implies that this star is 



more than 99% depleted, then our calculations imply that it must be older than 100 Myr for 
normal composition {X = 0.7), regardless of its effective temperature or mass (i.e., regardless of 
an evolutionary track.) Basri et al. (1996) reached the same conclusion by using the evolutionary 
tracks of Nelson et al. (1993). That our estimates agree is a consequence of the age constraint 



■^Similarly, a detection of lithium coupled with a lower bound on Tcff yields the maximum age. 
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Table 1. Radius and Mass of Lithium Depleted Pre-Main Sequence 

Stars 



Observed Inferred 



Star 


7<,,r(K) 


log(L/L ,) 


R/E. 




R(^f. 


HHJ 10 


3120 


-2.8 


0.14 


0.074 


1,2,(6) 


HHJ 36 


2825 


-2.5 


0.24 


0.11 


3,2,(6) 


LHS 248 


2775 


-2.5 


0.24 


0.11 


4,(1) 


GL 406 


2600 


-3.0 


0.16 


0.078 


5,(7) 


GL 569B 


2773 


-3.1 


0.12 


0.069 


5,(7) 



Note. — The 5% uncertainties in Teft correspond to a 10% uncertainties in radii and 
uncertainties of no greater than 10% in the minimum masses. The hthium non-detection 
is reported in the reference in parentheses. For stars with two references, the first is for 
effective temperature; the second, for himinosity. It is important to note that we have 
not calibrated the effective temperatures to a common scale. To quote our anonymous 
referee, "the pedigree of any temperature derived from observations cannot be ignored." 
The systematic errors between different Tgs scales can be as much as 10%. 

References. — (1) Martin et al. 1994b; (2) Stauffer et al. 1995; (3) Steele et 
al. 1995; (4) Besscl & Stringfellow 1993; (5) Burrows, Hubbard, & Lunine 1989; (6) 
Oppenheimer et al. 1996; (7) Magazzii et al. 1993 
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(in this low-mass range) depending mostly on L and being nearly independent of Tefj. As noted 
by Basri et al. (1996), this age agrees with the age of the Pleiades inferred from main sequence 
turn-off with convective overshoot (Mazzei & Pigatto 1989; Maeder &: Meynet 1991). In contrast, 
to share the 60-70 Myr age of the Pleiades inferred from main sequence turn-off without convective 
overshoot ( Mazzei &: Pigatto 1989| ), HHJ 3 must have either a hydrogen content X ~ 0.5 (very 
unlikely) or a luminosity 50% higher than that reported. 
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Fig. 2. — An example of using the formulae for depletion time (eq. radius (eq. y]), and luminosity (eq. Q) 

to bound the age of the Pleiades. For both HlfJ 3 and Calar 3, we plot a contour of constant lithium abundance 
(dashed for HHJ 3; dotted for Calar 3) and two constant luminosity contours (solid lines). These constr aints hound 
for e ach star a region of M-t space (shaded areas '). HHJ 3 has an effective temperature Teff = 2800 K ( Steele et al 
1995| ), a luminosity -2.83 < log{L/LQ) < -2.73 ( Basri et al. 1996[ ), and [Li] < 1.16. For Calar 3, the region of M-t 
space is bounded by Teff = 2600 K, —3.21 < log(L/-LQj < — 3.U1, and [Li] > 2.5. The overlap in t of the shaded areas 
sets a minimum age (105 Myr) and maximum age (not shown) of the Pleiades. Within this age range, each star then 
has a minimum and maximum allowed mass. (Rebolo et al. 1996). 



If both L and Teg are known, the absence (presence) of lithium sets a meaningful lower (upper) 
limit on both the mass and the age of an individual star. This method had its first practical 
application in the work of Basri et al. (1996), who used it to find a new age for the Pleiades and 
argued that the Pleiades' faintest members are substellar objects. Instead of plotting isochrones 
and evolutionary tracks on an L-Teff diagram, we find it more convenient to directly relate the 
two unknown quantities, mass and age, on an M-t plot (Figure 2). On this plot, the observed 
luminosity is represented by a line with positive slope found from the contraction relation (Q), 
while the limit on lithium abundance is represented by a line with negative slope inferred from the 
depletion formula (pA[). For a given Tgfj, the star must lie in a swath of the M-t plane defined by 
the estimated luminosity range. The lithium abundance constraint intersects this swath, setting 
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limits on age and mass. 

As an example, the hatching in Figure ^ denotes the allowed regions of the Pleiads HHJ 3 and 
Calar 3 (which shows lithium and appears undepleted) (Basri, Marcy, & Graham 1996; Rebolo et al. 
1996). Assuming that the stars are coeval, Calar 3 must have a mass < 0.075Mq, and HHJ 3 must 
have a mass > O.OSMq, if the effective temperatures ( pteele et al. 1995 ) are correctly measured. 



The recent detection of lithium in the Pleiades star PPl 15 (log[A'^Li/-^H] ^ —10.84 [ pasri et al. 



19961 ), which has T^s = 2800 ± 150 K and log[L/L0] = -2.80 ± 0.10 ( |Rebolo et al. 19961 ), provides 



a better upper bound on the age of the Pleiades than Calar 3. In particular, we find that the 
maximum age of PPl 15 (and hence of the Pleiades) is 145 Myr while its mass is constrained to 
lie in the range O.O7M0-O.O9M0. 



5. Conclusions 

Assuming that a contracting star is fully mixed and that the time it spends prior to Hayashi 
contraction is negligible compared to typical depletion times, we have derived simple analytical 
relations (equations @, @, |ll|, and |12|) for the radius, luminosity, and age of a star as a function 



of lithium depletion. These formulae demonstrate the dependence of the time of lithium depletion 
on the mass, composition, and effective temperature of the star. Reasonable agreement with other 
published theoretical calculations supports the use of our results to evaluate observations. 

We outline a method for using the observed T^g and L of both depleted and undepleted stars 
in a cluster to constrain its age. This method (complementary to that used by Basri et al. [1996] 
to date the Pleiades) is also applicable to clusters of age 10-100 Myr, where higher mass stars 
(up to O.5M0), are presently depleting lithium. These stars are relatively luminous (L ^ 10~'^Lq) 
when depleting lithium and are hence more easily observed. For example, a O.3M0 star will have 
depleted 50% of its lithium at about 16 Myr, when it has a luminosity O.O4L0 (for Tefj = 3300 K). 
Because stars with a range of masses will deplete simultaneously (the mass of a 99% depleted star 
is approximately 1.4 times the mass of a 50% depleted star), observations of stars within this mass 
range will constrain the cluster's age and the relative ordering of stellar masses. In a young cluster, 
the most stringent bounds are obtained by observing the dimmest depleted star and the brightest 
undepleted star. 

We have assumed that a fully convective star mixes material from the core to the photosphere 
faster than it contracts. Until a fully convective pre-main sequence star is observed that has both 
lithium and a central temperature hotter than the maximum burning temperature (4.4 x 10^ K for 
a O.5M0 star), this approximation goes unchallenged. In this sense, we have calculated the earliest 
time for depletion in the contracting star. Incomplete mixing can only delay the time of depletion 
to a given level. 
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